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I. INTRODUCTION 



Neutrinos have mass and non-trivial mixing angles, as neutrino oscillation experiments 
have shown . If neutrinos are Majorana particles, their mass at low energy is 

described by a unique dimension-5 operator Q| 

m u = L { HL){HL). (1) 



Using on 
operator 



y renormalizable interactions, there are only three tree-level realizations of this 



7]. The first one is the exchange of a heavy fermionic singlet. This is the celebrated 
seesaw mechanism jg[ 0, fiol, w hich we will call seesaw type-I. The second possibility is the 

as seesaw type-II. And lastly, 

one could also add one (or more) fermionic triplets to the field content of the SM uM- This 

n 

is called seesaw type-Ill in [7j, although this nomenclature is not universally accepted . 
The dimension-5 operator of eq. (CEJ) could also be generated at loop level. As the classical 



examples for loop generated neutrino masses we only mention the Zee model [17J (1-loop) 
and the Babu-Zee model 18] (2- loop), although many more models exist in the literature. 
A list of generic 1-loop realizations of eq. ([1]) can also be found in Q]. 

At "low" energies one can neither decide whether tree-level or loop physics generates eq. 
([I]), nor can any measurements of neutrino angles, phases or masses distinguish between 
the different tree-level realizations of the seesaw discussed above. Observables outside the 
neutrino sector are needed to ultimately learn about the origin of eq. ([ID. For loop generated 
neutrino masses, / in eq. ([1]) can be a very small number and the scale A at which new 
physics appears can be quite low, probably accessible at future accelerators such as the LHC 
or an ILC The "classical" tree-level realizations of the seesaw, unfortunately, can not be 
put to the test in such a direct way. This can be straightforwardly understood by inverting 
eq. (P), which results in A ~ f(^^) 10 15 GeV. 

Indirect inside into the high-energy world might be possible in supersymmetric versions of 
the seesaw. In the renormalization group equations for the soft SUSY breaking slepton mass 
parameters terms proportional to the neutrino Yukawa couplings appear. If the scale where 



1 Barr and Dorsner 



14| . for example, add additional singlets to the seesaw type-I. This version of the seesaw 



which the authors call type-Ill - might be named "double seesaw, variant-II" to distinguish it from the 



original double seesaw 15j, see also the related work in 16 1. 
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the right-handed neutrinos and/or the triplet decouple is below the scale at which SUSY 
breaks, lepton flavour violating (LFV) entries in the Yukawa matrices then induce LFV off- 
diagonals in the slepton mass matrices. This effect potentially leads to large values for lepton 
flavour violating lepton decays, such as \i — > e + 7, even if the soft masses are c omp letely 
flavour blind at high scale, as was first pointed out for the case of seesaw type-I in 19[ . It is 
maybe not surprising then that with the increasingly convincing experimental evidence for 
non-zero neutrino masses a number of articles have studied the prospects for observing LFV 
jrocesses, both at low en erg ies and at future colliders, within the supersymmetric seesaw 
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Despite the fact that a minimal seesaw type-II has fewer free parameters than the seesaw 
type-I, type-I seesaw has received considerably more attention in the literature. Prob ably 



one of the reasons for this preference is gauge coupling unification. As is well known [34, 35], 
the SM gauge couplings unify within the minimal supersymmetric standard model (MSSM) 
at a scale around Mq — 2 x 10 16 GeV, if the SUSY particles have masses around the electro- 
weak scale. Adding gauge singlets does not destroy this nice feature of the MSSM. However, 
a scalar triplet with mass below the GUT scale changes the running of g\ and <?2 in an 
unwanted way and gauge coupling unification is lost |36J. A simple way to cure this defect 
of the seesaw-II consists in adding only complete SU(5) multiplets (or GUT multiplets which 
can be decomposed into complete 577(5) multiplets) to the standard model particle content. 
In this way the scale where couplin gs u nify remains the same (at one loop level), only the 
value of the GUT coupling changes [371]. 

In this paper we calculate lepton flavour violating branching ratios of the scalar tau as well 
as LFV lepton decays at low energies, such as /j — > lj + 7 and Zj — > 31 j. For definiteness, we 
assume minimal Supergravity (mSugra) boundary conditions and fit the observed neutrino 
masses by a seesaw mechanism of type-II. We will discuss two different realizations. The 
first one is based on adding one pair of triplets to the MSSM, from which only one couples 
to the standard model leptons. This is the simplest supersymmetric version of the type-II 
seesaw. The second model we consider consists in adding a pair of 15 and 15 multiplets 



361 ] . This second option allows to maintain gauge coupling 



to the MSSM particle content 
unification also for M 15 <C Mq. 

We compare the sensitivities of low-energy and accelerator measurements and study their 
dependence on the unknown seesaw and SUSY parameters. Absolute values of LFV stau 
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decays and LFV lepton decays depend very differently on the unknown SUSY parameters. 
For a light SUSY spectrum, say slepton masses below 200 GeV, the current upper bound on 
Br(/i — > e + 7) limits seriously the possibility to observe LFV scalar tau decays. However, 
for heavier sparticles low energy data very rapidly looses its constraining power and large 
LFV at the LHC is allowed by current data. 

While absolute values of LFV observables depend very strongly on the soft SUSY breaking 
parameters, we discuss how ratios of LFV branching ratios can be used to eliminate most of 
the dependence on the unknown SUSY spectrum. I.e. ratios such as, for example, Br(f 2 — > 
e + X?)/Br(f2 — > fi + Xi) are constants for fixed neutrino parameters over large parts of 
the supersymmetric parameter space. Measurements of such ratios would allow to extract 
valuable information about the seesaw parameters: In the minimal type-II seesaw case these 
ratios can be calculated as function of measurable low-energy neutrino data. For the more 
involved case of the 15 + 15 model this simple connection is lost in general, but relations to 
neutrino data can be (re-) established in some simple, extreme cases for the Yukawa matrix 
Y15. We therefore study such ratios in some detail, first analytically then numerically. 

The presence of new non-singlet states below the GUT scale does not only affect the 
running of gauge couplings but also the evolution of the soft SUSY breaking parameters. 



Measurements of soft SUSY masses at the LHC and at a possib 



28, 



38 



LC therefore contain indi- 
|. ^From the different soft 



rect information about the physics at higher energy scales 
scalar and gaugino masses one can define certain "invariants", i.e. parameter combinations 
which are nearly constant over large ranges of the mSugra parameter space 39|, at least 
in leading order approximation. If the measured values of all the invariants depart from 
the mSugra expectation in a consistent way, one could gain some indirect estimate of the 
mass scale of the new particles, the scale of the seesaw type-II. We discuss first some lead- 
ing order analytical approximation, before showing by numerical calculation the limitations 
of the simplified analytical approach. While the different invariants indeed contain useful 
information about the high energy physics, reliable quantitative conclusions about the mass 
scale of the 15 require highly precise measurements of soft masses as well as a full numerical 
2-loop analysis. 

The rest of this paper is organized as follows. In the next section we will recall the 
basic features of the supersymmetric seesaw type-II and discuss a 577(5) motivated variant, 
which adds a pair of 15 and 15. Section (Till then discusses analytical solutions for the RGEs 
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and presents estimates for slepton mixing angles and the corresponding LFV observables. In 
Section [IV] we present our numerical results for LFV decays at low energies and accelerators. 
This numerical study demonstrates the reliability of our analytical approximations for the 
LFV observables. We then discuss soft masses and the seesaw type-II scale, demonstrating 
by a numerically exact calculation that for soft masses the leading order approximations are 
not accurate enough to draw quantitative conclusions. We then summarize in section |V] 



II. SETUP: MSUGRA WITH SEESAW TYPE II 



In this section, to set up the notation, we briefly recall the main features of the seesaw 



type 

of |36| 



I and mSugra. We then outline a simple S77(5) motivated model based on the work 



A. Supersymmetric seesaw with triplet (s) 

In supersymmetry at least two SU{2) triplet states Ti j2 with opposite hypercharge are 
needed to cancel anomalies. Thus, the minimal SUSY potential including triplets can be 
written as 

W = Wmssm + ^= (Y^LiTxLj + AiifiTWi + X 2 H 2 T 2 H 2 ^ + M T T{T 2 . (2) 

Here T\ [T 2 ] are supermultiplets with hypercharge Y = 1 (Y = — 1) and H\ 2 are the 
standard Higgs doublets with Y = =Fl/2. The matrix Yp is complex symmetric, Ai^ are 
arbitrary constants and Mt gives mass to the triplets, supposedly at a very high scale. Note 
that only T\ couples to the SM leptons, thus in the minimal (supersymmetric) model with 
two triplets the only source of lepton flavour violation resides in the matrix Yj*. 

Integrating out the heavy triplets at their mass scale the dimension-5 operator of eq. flTJ 
is generated and after electro-weak symmetry breaking the resulting neutrino mass matrix 
can be written as 



m,. 



v\ A 2 



Y T . (3) 



2 M T 

where v 2 is the vacuum expectation value of Higgs doublet H 2 and we use the convention 
(Hi) = Note that eq. (j3J) depends on the energy scale. m u is measured at low energies, 
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whereas for the calculation of m v we need to know A2, Yt and Mt as input paramters at the 
high scale. One can use an iterative procedure to find the high scale parameters from the 
low energy measured quantities, as explained in section IIVI In the basis where the charged 
lepton masses are diagonal, eq. (j3J) is diagonalized by 



m. 



U T -m v - U, 



(4) 



where the neutrino mixing matrix U is, in standard notation 
/ 

U = 



401 ]. given by 



C12C13 



S12C13 



V 



"S12C23 — Ci2S23Si3e l<5 C12C23 — Si2S23 s 13 e4<5 



si 3 e 

S23 c 13 



S \ ( giai/2 



\ 








(5) 





e ia 2 /2 o 

S12S23 - ci2C23Si 3 e JO -C12S23 - si2C23Si 3 e 10 C23C13 J y U 1 J 

Here Sjj = sin^ (cy = cos 6^). For Majorana neutrinos, U contains three phases: 5 
is the (Dirao) CP violating phase, which appears in neutrino oscillations, and 0:1,2 are 
Majorana phases, which can only be observed in lepton number violating processes. Neutrino 
oscillation experiments can be fitted with either a normal hierarchical spectrum (NH), or 
with inverted hierarchy (IH). If one does not insist in ordering the neutrino mass eigenstates 



m. 



1,2,3 with respect to increasing mass, the matrix U can describe both possibilities 



without re-ordering of angles. In this convention, which we will use in the following, m v , ~ 



[m 



V3 



0) corresponds to normal (inverse) hierarchy and S12, S13 and S23 are the solar (s©), 



reactor (sr) and atmosperic angle (sAtm) f° r both type of spectra. 
Note that 



U T -Y T -U 



(6) 



i.e. Yt is diagonalized by the same matrix as m v . If all neutrino eigenvalues, angles and 
phases were known, Yt would be fixed up to an overall constant which can be easily estimated 
to be 



M 7 



10 1& GeV 



0.05 eV 



At this points it might be worth recalling the main differences between seesaw type-II and 
seesaw type-I. In seesaw type-I there is one non-zero mass eigenstate for the light neutrinos 
for each right-handed neutrino added to the model. In contrast, seesaw-II can produce three 
non-zero neutrino masses with only one triplet. Thus the minimal model for seesaw type-II 
with only one triplet coupling to L has less parameters than seesaw type-I. We can count 
the new parameters in eq. (T5]): Yt being complex symmetric has 9 parameters. Additionally 
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we have Ai^ and Mt- All three could in principle be complex. However, field redefinitions 
on Ti and T2 can be applied to remove two of the three phases, thus there is a a total of 
13 parameters. Note, however, that only 11 of them are related to neutrino physics. Since 
we have the freedom to write down eq. fl2j) in the basis, where the charged lepton mass 
matrix is diagonal, we only have to add three charged lepton masses to the counting of free 
parameters. 2 This number should be compared to the 21 free parameters in seesaw type-I 



for three right-handed neutrinos 41]. At low energies a maximum of 12 parameters can 
be fixed by measuring lepton properties: 3 neutrino and 3 charged lepton masses, 3 angles 
and 3 phases. Thus from neutrino data neither seesaw type-II nor seesaw type-I can be 
completely reconstructed. However, especially important in the following is the fact, see eq. 
(JIJJ), that low-energy neutrino angles are directly related to the high-energy Yukawa matrix 
in seesaw-II, whereas no such simple connection exists in the seesaw type-I, see also the 
discussion in 22] . 



B. SU(5) inspired model with 15+15 

In this section we outline the basics of an SU(5) inspired model, which adds a pair of 15 
and 15 to the MSSM particle spectrum 36|. Our numerical calculations will all be based 
on this variant, since it allows to maintain gauge coupling unification for Mt -C Mq, as 
discussed in the introduction. 

Under SU(3) x SU L (2) x U(1) Y the 15 decomposes as 

15 = S + T + Z (8) 

S ~ (6,1,-^), r~ (1,3,1), Z~(3,2,i). 

3 

T has the same quantum numbers as the triplet T\ discussed above. The SU(5) invariant 
superpotential reads as 

W = -^Y 15 5 • 15 ■ 5 + ^Ai5h • 15 • 5* + -^\ 2 5h • 15 • 5 H + Y 5 10 -5-5* (9) 
a/2 a/2 a/2 

+ Y10IO • 10 ■ 5 H + M 15 15 • 15 + M 5 5 H • 5^ 



2 In the non-supersymmetric version of seesaw-II j^p — > j^t, with \i having dimension of mass, but the 
number of parameters related with neutrino physics does not change. 
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Here, 5 = (d c ,L), 10 = (u c ,e c ,Q), 5h = (tjH^) and 5# = (t,Hi). Below the GUT scale in 
the S77(5)-broken phase the potential contains the terms 

-^=(Y T LTiL + Y s d c Sd c ) + Y z d c ZL + Y d d c QH 1 + Y u u c QH 2 + Y e e c LH x (10) 
v2 

+ -^=(XiHiTiHi + \ 2 H 2 T 2 H 2 ) + M T T X T 2 + M^Z^ + M^ft + ^ X H 2 
v2 

The first term in eq. ({TO]) is responsible for the generation of the neutrino masses in the same 
way as discussed for the triplet-only case in the previous subsection. Y d , Y u and Y e generate 
quark and charged lepton masses in the usual manner. However, in adddition there are the 
matrices Yg and Yz, which, in principle, are not determined by any low-energy data. In 
the calculation of LFV observables in supersymmetry both matrices, Yj* and Yz, contribute. 
For the case of a complete 15, apart from threshold corrections, Yr = Y$ — Yz- One can 
recover the results for the simplest triplet-only model, as far as lepton flavour violation is 
concerned, by putting Y s = Y Z = 0. 

As long as Mz ~ Ms ~ Mt ~ M15 gauge coupling unification will be mantained. The 
equality need not be exact for successful unification. In our numerical studies we have taken 
into account the different running of these mass parameters but we decouple them all at the 
scale M T (M T ) because the differences are small. 



III. ANALYTICAL RESULTS 



A. Approximate solutions for the RGEs 



In mSugra one has in total five parameters at the GUT scale (42j. These are usually 
chosen to be Mo, the common scalar mass, M1/2, the gaugino mass parameter, Aq, the 
common trilinear parameter, tan j3 — ^ and the sign of /i. For the full set of RGEs for the 
15 + 15 see 36[. In the numerical calculation, presented in the next section, we solve the 
exact RGEs. However, the following approximative solutions are very helpful in gaining a 
qualitative understanding. 
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The gauge couplings are given as 



5a em (m z ) a em (m z ) 

Ui(jn z ) = — — , a 2 {m z ) = — — — , 11 

3 cos^ Vw sin 9 W 

oci(m z ) 

ai{rn SUSY ) = — 



i(rnz) h SM i ncr m susY 
®i{m SUSY ) 



Oi{M T ) = - 

1 - log J^- 

<Xi(M G ) = - 



SUSY 

Oi(Mr) 



^i + A^log- 



b SM = (b h b 2 , h) SM = (g, -f , -7) for SM and b = (b u b 2 , b 3 ) MSSM = (f , 1, -3) for MSSM. 
Mt denotes the mass of the triplet (15-plet). For the case of the complete 15-plet one finds 
Abi = 7 whereas for the case with triplets-only one finds Ab\ = 18/5, Ab 2 = 4 and A63 = 0. 
Using the equality ai(Ma) = a 2 (Mc) determines the GUT-scale Mq via 

lo M c 1 (12) 

° g Mf ai(m S usY)a 2 (m S usY)(b 1 + A61 - b 2 - Ab 2 ) 



2 

m SUSY 



M 2 \ 

4:TT(a 2 (m SU sY) ~ oi^itisusy)) + ai(msusY)a2(msusY)(b 2 - b{) log _ 2 T J 



Note, that in the case of the complete 15-plet Mq is independent of Mt- For the gaugino 
masses one finds 

MtimsusY) = °^m^ M l/2 . (13) 
a\Ma) 

Eq. ({TBI) implies that the ratio M 2 /M\, which is measured at low-energies, has the usual 
mSugra value, but the relationship to My 2 is changed. Neglecting the Yukawa couplings 
Y15, see below, for the soft mass parameters of the first two generations one obtains 

H< + M(^)V/Am;%, (14) 



m 

j 

1=1 



a x (M T ) M T 
-h log 



n 



h + A6. 

The various coefficients c{ are given in table [TJ 



Air m 
i a(M G ) n . , . , Ml 



T 



(15) 



9 
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E 


L 


D U Q 


J 


6 


3 


2 8 1 




5 


10 


15 15 30 


c 2 





3 
2 


| 


c 3 








8 8 8 
3 3 3 



TABLE I: Coefficients c{ for eq. (H 
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10 11 10" 10 1J 10" 10" 10" 
M n = M T [GeV] 

FIG. 1: Four different "invariant" combinations of soft masses (left) versus the mass of the 15-plet, 
Mi 5 = Mr. The plot assumes that the Yukawa couplings Y15 are negligibly small. The calculation 
is at 1-loop order in the leading-log approximation. 

Individual SUSY masses depend strongly on the initial values for M and M\i2- However, 
one can form different combinations, such as 



{m\-ml)/Ml 



a(M G ) V 



«2(m T ) V , 

— 7 r J 2 + J 2 

a{m G ) J 
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To 



ai{m T ) \ 
a{m G ) ) 1 



ai(m S usY)J 

which, to first approximation, are constants over large regions of mSugra space. We will call 
such combinations "invariants". 

Figure ([!]) shows four different invariants as a function of M i5 = M T , calculated using eqs 
(fT3]) - (|14|) . For Mt = Mq one reaches the mSugra limit. For lower values of Mt one obtains 
a logarithmic dependence on the value of Mt- If all the different invariants depart from their 
mSugra values in a consistent way, measurements of these parameter combinations can be 
used to obtain indirect information about the seesaw scale. In practice the "invariants" do 
depend on the SUSY spectrum and thus, indirectly still depend to some degree on the initial 
values of M and Mi/ 2 . We will discuss this point in more details in the numerical section. 
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For the off-diagonal elements of the slepton mass matrix, we will discuss only the left 



sector, since right slepton mass parameters do not run to first order approximation 



36|. In 



our numerical calculation we do solve the RGEs exactly and confirm this expectation. Off- 
diagonal elements are induced in m? due to the non-trivial flavour structure of the matrices 

Li 



Yt and Yz- Yt an d Yz appear symmetrically in the RGEs [36|]. Since only Yp can be fixed 
from low-energy data, for a general Yz the off-diagonal entries of m| do not follow any 
correlation with low-energy physics. For this reason in the following we will consider two 
extreme cases: (a) Yz = Yp, we will call this the 15-plet case; and (b) Yz = 0, we will refer 
to this as the triplet case. 

For m~ one finds the following approximation in the case of the 15-plet: 



lo M C } 

+3(A -^M[(t)- 7 -M>(t)f 

+3{Ao ~ 2^4 M{(t) ~ l m) ~ t M ' 3{t))2 ) dt (16) 

m) = M ^ i 1 ~ i + ^/a W mJ < 17 > 

In case of the triplet one finds 



! l/2 



+ 3 (A -^Mi(t)- 7 -Mi(t)Adt (18) 

The integration over t can be done analytically leading to corrections to the formulas for 
(Am|) y . 

We have found that the approximation formulas shown above work less well than the 
corresponding formulas for the seesaw type-I case and only give a rough order of magnitude 
estimate. The reason for this difference is, that in seesaw type-I the Y v hardly run, unless 
left neutrinos are very degenerate, as either the Yukawas themselfes are small or in the case 
of large Yukwas the contribution from gauge and (top) Yukawa couplings nearly cancel each 
other. Such a cancellation does not take place in case of Yt, thus leading to significantly 
stronger dependence of Yt on the renormalization scale and consequently larger differences 
between the numerical solutions and eqs (JT6|) and (TT8"}) . 
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However, we have found that it is possible to improve the accuracy of the approximation 
formulas using the results of the next subsection, see eq. fl25|) . The idea here is to replace 
the running Yukawa coupling Yp by the measured low-energy neutrino masses and angles 
times the unkonw coupling A2. In case A2 is sufficiently small, this parameter does run very 
little and eqs fll6p and (I18p agree already very well with the numerical results. For large 
values of A2, we can find approximate solutions for the RGE for this parameter, following 



28]. 



the procedure outlined in 
We define 

A 2 v 2 

X = ^, Y t = ¥-. (19) 

M 2 

The solution for the RGE for A2 is then given in terms of X(t) by (t = log -Jfr ) 

X(t) = X ( M G>x(t) m 



u x {t') 



l + lY t (Mg)t'Y^ 

Exit') 



1 9 



6l + Abi , „ r .A 5?* A-i-r^ / b2 + Ab 



/(61+A61) / 1 1 Ah \ ^7/(6 2 +A6 2 ) 



E x (t') = I 1+ ^ ^(Mg^J ( 1+ ^ a 2 {M G )t' 

we have found that, assuming an approximately constant Yj, the above equations become 
easy to solve and describe the running of A2 to a rather good approximation. Eq. fl20l and 
eq. (|25l) . together with eqs (Tl6|) and (TT8T) then allow to estimate LFV entries in m| up to 
an accuracy of typically a few percent. 



B. Analytical results for flavour violating processes 

Here we concentrate exclusively on the left-slepton sector. Taking into account the dis- 
cussion given above this is expected to be a reasonable first approximation. The left-slepton 
mass matrix is diagonalized by a matrix R l , which in general can be written as a product 
of three Euler rotations. However, if the mixing between the different flavour eigenstates 
is sufficiently small, the three differenct angles can be estimated by the following simple 
formula 

(Am|)y 

% " (Am|)« - (Am|),-/ (21) 
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LFV decays are directly proportional to the squares of these mixing angles as long as all 
angles are small. Taking the ratio of two decays, for example, 



5r(f 2 -+e + X ?) r^efV f( Am l 



Br{? 2 -+fi + Xi) V ^ ' V(Am 



113 



(22) 



one expects that (a) all the unknown SUSY mass parameters and (b) the denominators of 
eq. (l2Ti) cancel approximately. To calculate estimates for different ratios of branching ratios 
we define 

r Z = (23) 



| (Am 



\kl\ 



where the observable quantities are {r^) 2 . Of course, only two of the three possible combi- 
nations that can be formed are independent. 

We next derive some analytical formulas for (r^) 2 in terms of observable neutrino param- 
eters. The neutrino Yukawa coupling Y? can be written in terms of observable parameters 

Y T = ^^-m u = ^^-U* ■ diagfmi, m 2 , m 3 ) ■ W. (24) 
V2X2 V2X2 

The running of the soft-SUSY breaking slepton mass matrix is proportional to the 

parameter combination (Y^Yt)^. 3 This combination can again be expressed in terms of 
low-energy neutrino observables times an unknown scale: 



T)ij 



2M T l2 

2 



m 



(U ■ diag(m 2 , m 2 ., m 2 ) • U ] ).. 
J2 U ikU* k ml, (25) 

k 

The different off-diagonal entries are explicitly given as 

(Y}Y T ) 12 = m" 2 [UnU^ml + tA 2 £/ 22 m 2 + ^13^23^3] > (26) 
(Y}Y T ) 13 = m~ 2 [UnU^m 2 + U 12 U; 2 m 2 2 + U 13 U* 33 ml] , 
(Y±Y T ) 23 = m~ 2 [U 21 U* 3l ml + £/ 22 f/ 3 * 2 m 2 + U 23 U* 33 m 2 3 ] . 

Inserting the convention for the matrix U from eq. (0) results in 

(Y}Y T ) 12 oc Ci 2 si 2 ci 3 c 23 (m 2 - m 2 ) - Ci3Si 3 s 23 e"^{(m 2 - m 2 ) + c 2 2 (m 2 - m 2 )}, (27) 
(YjY T ) 13 oc ci 2 si 2 ci 3 s 23 (m 2 - m 2 ) - Ci 3 Si 3 c 23 e"* <5 {(m 3 - mj) + c^m 2 , - m 2 )}, 
{Y}Y T ) 23 oc s 23 c 23 ( (s 2 2 - c 2 2 )(m 2 - ml) + c 2 3 {(m 2 - m 2 ) + c 2 12 (m 2 2 - ml)} 



Si 2 ci 2 si 3 (c2 3 e 15 - sl 3 e l5 )(m 2 - m{). 



For the triplet-only case. For the 15 case we assume Yz — Yr at Mq, see the previous subsection. 
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Note, that the off-diagonals can be expressed as a function of mass squared differences only, 
i.e. there is no dependence on the overall neutrino mass scale. However, again note that 
eq. (127)) depends on the energy scale, see the discussion below eq. ((3]) and in section (lIVp . 
Also it is worth mentioning that with the convention of U from eq. (j5J) the Majorana phases 
cancel in eq. (127)) . 

As a starting approximation for the following discussion, let us assume that the lepton 
mixing matrix has exact tri-bimaximal (TBM) form 43] 



/ 



U = U' 



TBM 







1 

111 

VE V3 V2 
1 1 1 



(28) 



As is well-known, eq. (128]) is an excellent first-order approximation to the measured neutrino 



mixing angles 



44|. For these values eq. (1271) simplifies to 

\(Y}Y T )u\ oc ^Am 2 & 



(29) 



\(Y}Y T ) 13 \ cc ^Am 2 



0' 



\<y}Yi 



T)23 



oc -Am a 



Atirr 



The ratios (r\ 2 ) 



1 and (r 2 12 ) 



,13^ 
23 ) 



2 Amg 

3 Am 2 ,, 



result. 



For the general mixing matrix one can derive {r 1 ^) 2 using eq. (1271) . For the currently 
allowed ranges of the neutrino parameters, the most important unknown turns out to be 
S13, as fig. ([2]) demonstrates. In this figure (r^) 2 are shown as function of s 2 3 for tan 2 9a = 1 



and tan 2 6*0 = 1/2, as well as for the Am 2 fixed at their best fit point values 44] . Currently 
s 2 3 < 0.05 at 3 a c.l. rly strongly depend on the value of S13 and there exists a special value 
of S13, for which either (r^f) or (r^f) even vanish, due to a cancellation between the different 
terms in eq. (1271) . Note, however, that (r^f) and (r^f) can not vanish simultaneously. Note 



also, that for tan 2 #a 



and (rll) 2 are symmetric under the exchange of 5 = 



«-> 5 = 7r. Also, for non-zero values of S13 the results depend on the assumed hierarchy of 
the left neutrinos and the simultaneous exchange of the cases (normal hierarchy) NH IH 
(inverse hierarchy) and (5 — 0) (5 — tt) leads to the same values for the different rjy in 
case of tan 2 6 a = 1. 

Table (JTTJ) shows the currently allowed ranges for the (r^) 2 for S13 = and s™ x for differ- 
ent assumptions about the remaining neutrino parameters for the different cases of NH and 
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FIG. 2: Square ratios (rj|) (light blue line, full line), (^23) (blue line, dashed line) and (^23) ( re d 
line, dotted line) versus s 2 3 for NH (upper panels), IH (lower panels) for 5 = (left panels) and 
5 = tt (right panels). The other light neutrino parameters have been fixed to their b.f.p. values. 



Note, that for tan 9a = 1, (^23) anc ^ ( r 23) are symmetric under the exchange of 5 = <-> 5 = tt. 
Also the simultaneous exchange of NH <-> IH and (5 = 0) <-> (<5 = tt) leads to the same values for 



the different r^, in case of tan 2 



1. 



IH. These values serve to indicate the allowed variations for due to other parameters than 
S13. As stated above, the allowed variation on s 13 is most important for the "uncertainties" 



m r 



n \2 

kl) ■ 



However, also the current error bar on tan 2 8 a leads to a sizeable variation on r^. 



Since Am Atm and AttIq are now known to much better precision than the neutrino angles, 
their variation is much less important for the {r^) 2 , as table (JTTJ) demonstrates. 

Finally, recall that all results presented in this subsection are based on the assumption 
that one of the extreme cases, Yz = Yt or Yz = 0, is realized. The former corresponds 
to the SU(5) inspired model with a complete 15 of section fill Bp . whereas the latter cor- 
responds to the simplest triplet-only model discussed in section fill Al) . However, we stress 
that departures in the ratios of LFV branching ratios from the values calculated in this 
subsection should not be interpreted as "ruling out" the seesaw type-II. Rather they should 
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NH 


IH 


5 = 


5 = IT 


5 = 


5 = IT 


812 = 1/V3 

s 23 = l/y/2 

S13 = 


(^il) 2 


1 




1 






[2.8,7.4] x 10~ 4 




[2.8,7.2] x 10~ 4 




& 


[2.8,7.4] x 10~ 4 




[2.8,7.2] x 10~ 4 




sia = l/\/3 

S 23 = 1/V2 

si3 = 




[1.2,1.4] 


[0.71,0.81] 


[0.71,0.81] 


[1.2,1.4] 




[0.12,0.13] 


[0.091,0.096] 


[0.085, 0.093] 


[0.11,0.12] 


& 


[0.091,0.096] 


[0.12,0.13] 


[0.11,0.12] 


[0.085,0.093] 


S12 / 1/V3 

S 23 + 1/V2 
S13 = 


(r\l? 


[0.49, 1.94] 




[0.49, 1.94] 




(rllf 


[1.8, 12] x 10~ 4 




[1.8,12] x 10" 4 




(r\l? 


[1.8, 12] x 10~ 4 




[1.8,12] x 10" 4 




sia / I/a/3 

S23 / l/\/2 

513 = s^r 


(rg) a 


[0.63,3.0] 


[0.35, 1.7] 


[0.35, 1.7] 


[0.63,3.0] 




[0.094,0.18] 


[0.062,0.15] 


[0.060,0.15] 


[0.089,0.18] 


(^l) 2 


[0.061,0.15] 


[0.093,0.18] 


[0.088,0.17] 


[0.058,0.14] 



TABLE II: The parameters are given for several values of the neutrino mixing angles. Si3 ax 
is the experimentally allowed maximum value: (s^*) 2 = 0.050 at (3<r) c.l. NH and IH are 
normal and inverted hierarchy of neutrino masses, respectively. The intervals correspond to (3a) 
experimental allowed range of neutrino oscillation parameters: sf 2 = 0.26-0.40, s% 3 = 0.34-0.67, 
Am| = (7.1 - 8.3) x 10~ 5 eV 2 and Am 2 Atm = (2.0 - 2.8) x 10" 3 eV 2 . In the top two rows only the 
mass squared splittings are varied, while for the lower set also angles are allowed to vary. 

be interpreted in the sense that one has to go beyond minimal scenarios. 



IV. NUMERICAL RESULTS 



In this section we present our numerical calculations. All results presented below have 



been obtained with the lepton flavour violating version of the program package SPheno 45]. 
Calculations are done for the 15-plet case, using the assumption Y z = Y T at M G , as discussed 
above. Unless mentioned otherwise, we fit neutrino mass squared differences to their best fit 
values 44| and the angle to TBM values. Our numerical procedure is as follows. Inverting 



1(3 



the seesaw equation, see eq. one can get a first guess of the Yukawa couplings for any 
fixed values of the light neutrino masses (and angles) as a function of the corresponding 
triplet mass for any fixed value of A2. This first guess will not give the correct Yukawa 
couplings, since the neutrino masses and mixing angles are measured at low energy, whereas 
for the calculation of m u we need to insert the parameters at the high energy scale. However, 
we can use this first guess to run numerically the RGEs to obtain the exact neutrino masses 
and angles (at low energies) for these input parameters. The difference between the results 
obtained numerically and the input numbers can then be minimized in a simple iterative 
procedure until convergence is achieved. As long as neutrino Yukawas are not too close 
to one we reach convergence in a few steps. However, in seesaw type-II the Yukawas run 
stronger than in seesaw type-I, thus our initial guess can deviate sizeably from the exact 
Yukawas. Since neutrino data requires at least one neutrino mass to be larger than about 
0.05 eV, we do not find any solutions for > A2 • 10 15 GeV. 

We have implemented the effects of the additional triplets (15-plets) including the two- 
loop contributions to the RGEs for gauge couplings and gaugino masses, one-loop contri- 
butions to the remaining MSSM parameters and one-loop RGES for the new parameters in 
SPheno. For consistency we have also included 1-loop threshold corrections for gauge cou- 
plings and gaugino mass parameters at the scale corresponding to the mass of the triplet. 
The MSSM part is implemented at the 2-loop level and, thus, in principle one should also 
include the effect of the 15-plets consistently for all parameters at this level. However, the 
correct fit of the neutrino data require that either the triplet (15-plet) Yukawa couplings are 
small and/or that Mt is close to Mq implying that the ratio Mt/Mg is significantly smaller 
than Mc/rriz and thus one expects only small effects. 



A. Numerical results for LFV 



The analytical results presented in the previous section allow to estimate ratios of branch- 
ing ratios for LFV decays. For absolute values of the branching ratios, as well as for cross- 
checking the reliability of the analytical estimates, one has to resort to a numerical calcula- 
tion. Below we show results only for a few "standard" mSugra points, namely for SPS3 46] 
and SPSla' [47]. However, we have checked with a number of other points that our results 
for ratios of branching ratios are generally valid. 
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FIG. 3: Lepton flavour violating branching ratios versus Mt = -M15 for the standard mSugra 
point SPS3 for two values of A2. To the left A2 = 0.05, to the right A2 = 0.5. The plots show 
Br{li — ► lj + 7) (top) and Br(j2 — > e, [i + x?) (bottom). Ratios of the different branching ratios 
follow closely the analytical expectations. The regios excluded by the current upper limit on 
Br{^L — ► e + 7) is shown also in the lower plot. 

Fig. shows examples of LFV decays for the mSugra point SPS3 as a function of 
Mt = M15 for two different values of A2. The upper plots show Br(li — > /j + 7), while the 
lower ones show Br(f 2 — > e, [/, + xD- We have also calculated Br(li — > 3lj), but these are 



not shown in the plots, because they follow very well the approximate relation [24. 



48| 



Br(Zi -> L + 7) 3tt 



mr 4 / 



(30) 



All LFV branching ratios show a very strong dependence on the value of Mt and due to the 
stronger running of parameters in the seesaw type-II case, compared to the seesaw type-I, 
the dependence on the seesaw scale is stronger than in seesaw-I 49] . See also the discussion 
in section (JTTTJ) . 

For the calculation shown in fig. ([3]), we have fitted the neutrino angles to exact tri- 
bimaximal values. One sees that, as long as the different LFV branching ratios are small, 
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ratios of branching ratios are constants, which follow very well the analytical expectations. 
Currently the most important phenomenological constraints comes from the upper limit on 
Br{n -> e+ 7), Br(fi -»■ e + 7) < 1.2 ■ 10~ n Q. Note that the "dip" in Br(fi -> e + 7) is 
due to a level-crossing of selectron and smuon mass eigenstates. 4 For SPS3 one finds that 
this limit rules out Br{r 2 — > fi + Xi) larger than a few percent, the exact number depending 
on the unknown parameter A2. Fig. ([3]) to the left (right) shows results for A2 = 0.05 
(A2 = 0.5). Recall that neutrino physics fixes only However, note also that the 

upper limit on Br(f 2 — > fi + Xi) depends only weakly on A2. 
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FIG. 4: As rig. ([3]), but for the mSugra standard point SPSla' and for A2 = 0.5. For a slepton 
spectrum as light as expected for SPSla' Br(n — ► e + 7) rules out the possibility to observe large 
lepton flavour violating slepton decays at the LHC. 



It is well-known that absolute values of LFV branchin g ra tios depend very strongly on the 



SUSY spectrum, for example Br(fi — > e+7) oc l/m 8 SUSY (2l|. Since both left-sleptons as well 
as (lightest) neutralino and chargino are approximately a factor of two heavier for SPS3 than 
for SPSla', one expects that Br(fi — > e + 7) gives a strong constraint on the observability 
of LFV at the LHC for SPSla'. This is confirmed numerically, as shown in fig. (J2J), which 
shows Br{li — > lj + 7) and Br{f 2 — > e, \i + xT) as function of Mr = Mi 5 for the example 
of A 2 = 0.5. Given the current limit on Br(/i — > e + 7) one expects Br(f 2 — > fi + Xi) ~ 
(few) 10~ 4 . Note that again we have fitted neutrino angles to tri-bimaximal values in this 



4 In mSugra the left-selectron is usually slightly lighter than the left-smuon. In the mSugra plus seesaw case, 
for both typc-I and typc-II seesaw, the additional Yukawas change the running of the slepton masses. For 
large Yukawas (i.e. large Mt) fitting current neutrino data requires couplings such that the smuon mass 
runs faster to smaller values than the selectron mass. However, the splitting between seletron and smuon 
mass eigenstates is expected to be too small to be measurable in most parts of the parameter space, see 
the discussion in the next subsection. 
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calculation and that ratios of LVF branching ratios follow closely the analytical expressions. 



B. Sparticles Masses and seesaw scale 

As disdussed in the analytic section, the running of soft parameters allows, in principle, 
an indirect determination of the seesaw scale. In this section we discuss numerical results 
for the running of the "invariants" defined above. Although below we show plots only for 
the combination (m| — m 2 ^)/M 2 we have checked numerically that all invariants shown in 
fig. (pQ) behave qualitatively in the same way. 




10 11 10 12 10 13 10 14 10 15 10" 

M 1S = M T [GeV] 

FIG. 5: "Invariant" (m| — rrv^/M 2 , calculated with negligibly small Yukawa couplings for two 
mSugra standard points. The figure shows a comparsion of different calculation. The curve labeled 
"Analytic" uses the formulas presented in the previous section. 1-loop and 2-loop stand for exactly 
solved numerical calculations using 1-loop and 2-loop RGEs. Note the significant shift when going 
from 1-loop order to 2-loop order. 



Fig. (J5J) shows (m| — m 2 ^)/M 2 as a function of = M 15 for SPSla' and SPS3 comparing 
different calculations. This plot assumes that the Yukawas of the 15-plet are negligibly 
small, i.e. neutrino mass are not correctly fitted in this calculation. The black line is the 
analytical calculation based on 1-loop RGEs and the leading-log approximation with an 
assumed msusY — 1 TeV. The dotted lines are the numerically exact results for this invariant 
using 1-loop RGEs, while the full lines are the exact results using 2-loop RGEs. Obviously 
the "invariant" does depend to a certain degree on the mSugra point, as already pointed out 



in section fillip . However, we also find a considerable upward shift of (m| — m|,)/Mf , when 
going from the 1-loop to the 2-loop calculation. Since the dependence of (m| — m 2 ^)/M 2 
on the value of Mt is only logarithmic, even such a moderate change in the invariant is 
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important, if one wants to extract an indirect estimate on Mt from such a measurement. 
Note that for the point SPSla' the calculation stops at M15 ~ 10 11 ' 6 GeV, the lowest value 
of M 15 for which correct electro-weak symmetry breaking occurs. 

We have checked by an exact numerical calculation that the other invariants shown in 
section (lllip suffer from similar changes when going from 1-loop order to 2-loop. In other 
words, if one wants to learn about the seesaw scale from measurements of the soft masses, 
a careful analysis at 2-loop order will be necessary. Also note that, due to the logarithmic 
dependence on Mt, highly precise measurements will be necessary, especially if Mt is large, 
say M T > 10 12 " 13 GeV. 
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FIG. 6: "Invariant" (m~ — m|,)/M 2 calculated with Yukawa couplings fitted to neutrino data, for 



an arbitrary choice of A2 = 0.5. The calculation uses 2-loop RGEs. Results are shown for SPSla' 
and SPS3. Neutrino angles are assumed to have exact TBM values. 

Fig. shows (m| — m^)/M 2 calculated with Yukawa couplings fitted to neutrino data, 
for an arbitrary choice of A2 = 0.5. The calculation uses 2-loop RGEs and results are shown 
again for the mSugra standard points SPSla' and SPS3. For Mt low, say Mt < 10 13 GeV 
or so in this example, Yukawa couplings which explain current neutrino data are too small 
to induce any significant effect in the determination of (m| — m^)/M 2 . 

However, for larger values of Mt sizeable differences between fig. and fig. ([HD show 
up. First of all, for negligibly small Yukawas the calculation can vary Mt freely up to the 
GUT scale. If instead we insist to fit neutrino masses, such large values for Mt are not 
allowed. The downward turn in (m| — m 2 ^)/Mf is due to Yukawas, which if larger than 
0(0.1) contribute sizeable in the running of the soft parameters. In the example shown 
in this figure A2 = 0.5 has been chosen. For smaller values of A2 again for fixed values of 
the Yukawa couplings to fit neutrino masses a lower Mt is required. Correspondingly, for 
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smaller A 2 the effect of the Yukawas is seen for smaller values of Mt- 

It is also found that slepton mass parameters of the first and second generation run 
differently for large values of Mt, see fig. This difference can be traced to the fact that 
we have fitted neutrino angles to take exact TBM values. In this limit, m| oc Amg, while 
m| oc Am\ tm . Thus, at the largest values of Mt sizeable mass differences between 1st 
and 2nd generation sleptons show up. This difference is expected to be smaller for non-zero 
values of S13. Note, however, that for the example points shown in fig. (JSJ), there is the upper 
limit on Mt from Br(fi — > e+7), discussed in the last subsection. For SPSla' Mt < 1.5- 10 13 
GeV, for SPS3 Mt ^ 6 • 10 13 GeV. This limits the range of Mt where differences between 
1st and 2nd generation slepton masses might be observable. We mention that a recent 



paper 



50| claims that mass differences between smuons and selectrons can be measured 



very accurately, even at the LHC Depending on the mSugra point (m~ — m|)/(m~ + mf) 



as small as O(10 4 ) might be measurable 50j provided the leptons have sufficient energy to 
pass the experimental cuts. 
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FIG. 7: Branching ratios for LFV decays versus (m| — m|,)/M^ for SPS3 for two different values 
of A2. Measuring both types of observables allow in principle to disentangle A2 and Mt- 



All observables discussed so far are sensitive only to a combination of Mt and A2. If, 
however, both LFV decays as well as (m| — m 2 ^)/Mf could be measured in the future, one 
could disentangle the two parameters, in principle, by combining both measurements. This 
is demonstrated in fig. ([7]), which shows LFV decays, Br(fi — > e + 7) and Br{j2 —>■ e, fi + Xi) 
versus (m| — m 2 ^)/M 2 , for two different values of A2. Note again that the "dip" in Br(fi — > 
e + 7) is due to a level-crossing of selectron and smuon mass eigenstates. However, again 
we warn that a full 2-loop calculation is needed, before any quantitative conclusions could 
be drawn from such a measurement. 
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V. CONCLUSIONS 



We have studied phenomenological implications of the supersymmetric version of the 
type-II seesaw within mSugra. We have calculated lepton flavour violating observables, such 
as Br(li — > lj + 7) and LFV scalar tau decays. We have found branching ratios for LFV 
violating stau decays are large enough to be detectable at the LHC in principle. We have 
pointed out that in the simplest case of only one triplet coupling to the SM leptons, ratios 
of LFV branching ratios can be calculated from low energy neutrino data only. However, 
for the case of a complete 15 multiplet the situation is not as straightforward. In the 
SU(5) inspired model the Yukawa couplings Y T and Y z are related to the Y 15 and the 
conclusions remain unchanged. However, allowing and Yz to be free parameters, the 
relation with neutrino physics is lost. Thus, seesaw type-II can not be ruled out by any 
LFV measurements in general. Instead measuring ratios of LFV branching ratios can be 
understood as a consistency check for the minimal seesaw type-II models. 

We have also calculated the soft masses as a function of the seesaw parameters. As 
discussed in some detail, there are certain combinations of soft masses, which are approx- 
imately constants over large regions of mSugra space. These "invariants" contain indirect 
information about the seesaw scale. Measuring SUSY masses as precisely as possible will 
therefore allow to constrain the scale of seesaw type-II indirectly. However, theoretically 
there are many possibilities, why any single of the "invariants" we discussed could depart 
from the simplest mSugra expectations. Only a consistent departure of several "invariants" , 
together with measurements of LFV processes, could therefore be taken as a hint for seesaw 
type-II. 
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APPENDIX A: CONTRIBUTIONS TO THE (3 FUNCTIONS 



Using general formulas by 5l| we obtain for the RGEs of the gauge couplings: 
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